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<N : ABSTRACT 

c ■ 

^ • We compute the kinematical dependence of the peak luminosity, the pulse 

■ width and the spectral lag of the peak luminosity on the viewing angle 9^ of a 
^ ' jet. For appropriate model parameters we obtain the peak luminosity-spectral 

^ ■ lag relation similar to the observed one including GRB980425. A bright (dim) 

^ ■ peak with short (long) spectral lag corresponds to a jet with small (large) 

' viewing angle. This suggests that the viewing angle of the jet might cause 

^ ■ various relations in GRBs such as the peak luminosity-variability relation and 

^ ■ the luminosity-width relation. Our model also suggests that X-ray rich GRBs 

■ (or X-ray flushes or Fast X-ray transients) are typical GRBs observed from large 
p i' 9y with large spectral lag and low variability. 

o ■ 

^ ■ Subject headings: gamma rays: bursts — gamma rays: theory 

V , 

1. INTRODUCTION 

Afterglows of the GRBs are beheved to be produced by the external shocks while the 
GRB itself is by the internal shocks. There is a growing evidence that some of afterglows, 
such as GRB000926, GRB000301C, GRB990510 and GRB990123, are collimated with an 
opening half-angle of AO^""^ ~ 0.1. However, the collimation half-angle of the internal 
shocks remains unknown since we can observe only the angular region inside the relativistic 
beaming half-angle ~ where 7 is the Lorentz factor of the flow and is larger than ~ 100 
to avoid the compactness problem (e.g., Piran 1999). Since the angular size of a causally 
connected region is also ~ 7"^, the minimum possible A9 is ~ 7"^ ^ 10"^ ^ A9^"'\ that is, 
can be much smaller than A^'^'^^ 

We here assume that the internal shocks consist of such sub-jets with the collimation 
half-angle of A^ <C A9^"'\ If we happen to observe a bright sub-jet, the isotropic luminosity 
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becomes very large so that the sub-jets model can interpret (1) no correlation between 
the gamma-ray fluence and the redshift, (2) the weak correlation between the gamma-ray 
fluence and the afterglow flux as well as (3) the extremely energetic bursts ~ 10^^ ergs 
(Kumar & Piran 2000; but see also Freedman & Waxman 1999). The X-ray pre/post-cursor 
(Murakami et al. 1991; Sazonov et al. 1998; in't Zand et al. 1999) can also be explained in 
the sub-jets model, if some of off-axis sub-jets are ejected earlier (for precursor) or later (for 
postcursor) than the main sub-jets (Nakamura 2000). X-ray rich GRBs (Strohmayer et al. 
1998; Prontera et al. 2000; Heise et al. 2001; Kippen et al. 2001) are also expected in the 
sub-jets model if none of sub-jets point to our hne-of-sight (Nakamura 2000). The "sub-jet" 
is sometimes termed "cannonball" (Dar & De Riijula 2000a, b; Plaga 2000), "mini-jet" 
(Shaviv & Dar 1995) and so on. 

On the other hand, there are some possible relations between the peak luminosity, the 
spectral lag, the pulse width and the variability of GRBs. Prom six GRBs with known 
redshifts, Norris, Marani & Bonnell (2000) found that the isotropic peak luminosity is 
inversely proportional to the spectral lag which is deflned as the time lag of the peak 
luminosity between different energy bands. Penimore & Ramirez-Ruiz (2000) found that 
the luminosities of seven bursts with known redshifts are correlated with the variabilities 
of the bursts (see also Reichart et al. 2000). If both relations are correct, there should be 
a relation between the variability and the spectral lag. Recently Schaefer, Deng & Band 
(2001) showed that this is the case by plotting the variability and the spectral lag of 112 
GRBs. Therefore at present there is a strong support that the luminosity-spectral lag as 
well as the luminosity- variability relation are correct. 

There are several attempts to explain the luminosity-spectral lag and the luminosity- 
variabihty relations. Salmonson (2000, 2001) proposed that it is due to the difference of 
the observed cooling or deceleration time which depends on the variation in line-of-sight 
velocity among bursts. Schaefer (2001) interpreted the luminosity-spectral lag and 
luminosity- variability relation as a consequence of the conservation of energy when radiative 
cooling dominates and the dependence of the luminosity and variability on the Lorentz 
factor of the ejecta, respectively. However Wu & Penimore (2000) pointed out that the 
synchrotron cooling time scale is much shorter than the lag time scale. Plaga (2000) 
argued that the luminosity- variability relation might be obtained in the framework of the 
cannonball model for the given power low power density spectrum of the GRBs. 

In this Letter, under the assumption of short cooling time scale, we will show that 
depending on the viewing angle of a single jet to our line-of-sight, the luminosity-spectral 
lag and the luminosity-variability relation may arise kinematically. In the sub-jets model, 
we at present do not know the number of sub-jets in A^^"). If it is small, that is, the filling 
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factor of the sub-jets is small, we should study effects of many independent sub-jets based 
on the present results of a single jet. If the number of sub-jets is large, we may apply the 
present results regarding ~ A^^'^^ In any case it is important to study the viewing 
angle dependence of various physical quantities, which is the main purpose of this Letter 
irrespective of the details of the sub-jets model. 



2. INSTANTANEOUS EMISSION FROM AN EXPANDING THIN JET 

The internal shock occurs in the relativistic wind when the fast moving flow catches 
up the slow one. The wind can be modeled by a succession of relativistic shells (Kobayashi, 
Piran & Sari 1997; Daigne & Mochkovitch 1998). A collision of two shells produces a 
single pulse light curve, whose superposition makes a whole light curve. Here we consider 
a single pulse. There are three time scales that determine the temporal pulse structure: 
the hydrodynamic time scale, the cooling time scale, and the angular spreading time scale 
(Kobayashi, Piran & Sari 1997; Katz 1997; Fenimore, Madras & Sergei 1996). The coohng 
time scale is usually much shorter than the hydrodynamic time scale in the internal shocks 
(Sari, Narayan & Piran 1996), so that we may assume a delta function time dependence 
of the emissivity in the comoving frame. Let / and L be the width of the shell and the 
separation of the shells. Then the hydrodynamic time scale ~ //c and the angular spreading 
time scale ~ L/c determine the rise and the decay time of a single pulse, respectively. Since 
most observed pulses rise more quickly than they decay (Norris et al. 1996), for simplicity 
we may adopt an infinitesimally thin shell approximation as / ^ L. 

Let us consider an emitting thin shell that is confined to a cone of an opening half-angle 
A^^ and moves radially outward with the Lorentz factor of 7 = 1/ \/l — (3'^. A general formula 
to calculate the observed flux from an optically thin material is derived by Granot, Piran 
& Sari (1999) and Woods & Loeb (1999). Here we adopt their formulations and notations. 
Let us use a spherical coordinate system r — (r, 9, 0) in the lab-frame, where the ^ = 
axis points toward the detector and the central engine is located at the origin. Let also the 
detector be at a distance D from the source and a :— rsinO/D — ry/l — jjfl/D be the angle 
that a given ray makes with the ^ = axis. Then the observed flux at the observed time T, 

measured in erg s'^ cm'^ Hz'S is given by F,(T) = ^ ^^'^ d<f> J,^- aHa C.^^^] 
and = (1 — v' I'-^u)/ 13, where Q^, t — T -\- (rfi/c), am, and j^, are the direction towards the 
detector measured in the comoving frame, the lab- frame time, the maximum value of a, 
and the comoving frame emissivity in units of ergs s~^ cm~^ Hz~^ sr~^. Note here that a 
prime means the physical quantities in the comoving frame. 



If the emission is isotropic in the comoving frame, the emissivity has a functional form 
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of j'An'„ r, t) = AomSit - to)5ir - t,)H{M - \e - e,\)H [cos <t> - ( '-"^tU^inT' ^ 
where /(z^') represents the spectral shape and 9^ is the angle that the axis of the emission 
cone makes with the ^ = axis. The delta functions describe an instantaneous emission 
at t = to and r = ro, and H{x) is the Heaviside step function which describes that the 
emission is inside the cone. Then, the flux of a single pulse is given by 



Fu{T) 



2cAoroi'' A(l){T)f [1/7(1 - /3cos^(T))] 



L»2 



[7^(1-/3 cos e(r))r 



where 1 - /3cos^(T) = {c(3/ro){T - Tq) and Tq = to - r^jcii. For > and < d[T) < 
A9 - e^, A(f){T) = TT, otherwise A0(r) = cos-^ [^^^^^^SSw^] • < ^(^) ^^"^^ 

from to + while from 9^ - A9 to 9^ + A9 for 9j> A9. In the latter case, A0(r) = 
for 9{T) ^9„- A9. 



The observed spectrum of GRBs is well approximated by the Band spectrum (Band et 
al. 1993). In order to have a spectral shape similar to the Band spectrum, we adopt the 
following form of the spectrum in the comoving frame. 



I+ob 



1+ 




(2) 



where as {Pb) is the low (high) energy power law index, and s describes the smoothness of 
the transition between the high and low energy, as ^ —1 and Pb ~ —2.5 are typical values 
(Preece et al. 2000). Equations (1) and (2) are the basic equations to calculate the flux of 
a single pulse, which depends on nine parameters for 9^ <^ 1 and A^ -C 1: 71^0) l^v, 1^9, 
^o/c/37^ ^0 = ^0 - ro/cP, aB, Pb, s and 2cAoroj^/D^. 



3. PULSE PROFILE AND SPECTRUM 

In order to study the viewing angle dependence of a pulse from a single jet, we first 
fix other parameters: 7A^^ = 1, aB = —I, Pb = —2.5 and s = 1. The lab- frame frequency, 
the observed time and the flux per frequency are measured in units of 7z/q, ro/cP'-f'^ and 
2cAoro7^/-D^, respectively. The total fluence is given by S,^ = J^^"J^^ Fi,{T)dT, where 
Tstart = To + {ro/cp){l - p cos(max[0, 9^ - A9])) and T.^a = To + {ro/cP){l - P cos(^, + A^)) 
In Figure 1, we plot uSu as a function of the observed frequency 1/ by varying the viewing 
angle 7^^. As 7^^ increases, both the maximum frequency and uSu decrease. Here the 
maximum frequency lymax means the frequency at which most of the radiation energy is 
emitted. In Figure 1 we also show the analytically estimated dependence, Umax and 
{i^Si,)max oc 5~^, by triangles, where 5 — 7(1 — /9cos^„) ~ (1 + ^^7^)727 is the Doppler 
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factor (see loka & Nakamura 2001 for analytical estimates.). Figure 1 shows that if we 

observe the jet from O,, ^ several x 7"^ the maximum frequency Vmax is in the X-ray band 
while for 9^ ~ 7^^ it is in the gamma-ray region so that depending on the viewing angle, 
GRBs may be observed as X-ray rich GRBs in our model (Nakamura 2000). Note that if 
the emitting sub-jets distribute sparsely within the opening angle of the afterglow A^^"'' or 
if we are seeing the edge of the afterglow cone, the afterglow can occur as usual. 

Several pulse profiles are shown in the upper panel of Figure 2, where the observed 
frequency is taken as v = 200kcV(7t'Q/10'^kcV). The pulse profiles are similar to the 
FRED (Fast Rise Exponential Decay) shape observed in many bursts. Note here that 
the vertical axis is a log scale so that the straight line means the exponential decay. In 
the lower panels of Figure 2 we show the pulse width Wfwhm at half maximum and the 
product of the peak flux and the width i'F'P'^"-^ Wfwhm as a function of (1 -|- 7^^^) oc 5 
with the analytically estimated dependence on (1 -|- 7^^^) (see loka & Nakamura 2001). 
Figure 2 shows that the luminosity-width relation is given by z/F^*^"*^ oc Wp^y^j^ where 

+ 3 for 6^ ~ M and 7^^ ^ [10(7i/^/103keV)(z/„6,/200kcV)-^ - l]^/^^ r^^iis is 

consistent with the observation that the narrower pulse tends to be brighter in each burst 
following hiF'P^"-'^ oc Wp'^^j^^ (Ramirez-Ruiz & Fcnimore 2000). Our calculations predict 
that the exponent n ranges from k = — 2 -|- ~ — 3 to k = — 3 + 213b ~ —8 for dimmer 
bursts including X-ray rich GRBs which have the larger viewing angle in our model. 

4. PEAK LUMINOSITY-SPECTRAL LAG/ VARIABILITY RELATION 

In the upper panel of Figure 3, varying the viewing angle 7^^, we show the isotropic 
peak luminosity at frequency = 200keV(7z/Q/10'^keV) as a function of the spectral lag 
in an arbitrary vertical scale of 2cAoro'y'^ / D"^ . For simplicity, we define the spectral lag 
AT as the difference of the peak time between frequencies — 200keV(7z/o/10^keV) and 
ux — 20keV(7z/Q/10^keV), although the spectral lag is determined by the cross-correlation 
analysis in Norris, Marani & Bonnell (2000). Since there is no lag when the observer is inside 
the emission cone 9v < A^ in our present definition, all points with AT < 10~^(ro/c/97^) 
are plotted at AT = 10~^(ro/c/97^). In the lower panel of Figure 3 the solid line shows the 
corresponding viewing angle 7^^ as a function of AT. We adopt Tq/cP^'^ = 10 s. The peak 
luminosity first decreases with increasing spectral lag and turns back for large spectral lag 
since the different 9^ give the same lag. The analytical estimates, uF^'^'^'^ oc AT~^/^ and 
7^^ oc AT^/"^, shown by dotted lines agree well with the solid lines for 1 ^ 76*^ ^ 3 (see 
loka & Nakamura 2001). The deviation from the power law at low luminosity is due to 
the fact that the maximum frequency i^max becomes lower than the observed gamma-ray 
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frequency (see Figure 1). The observed luminosity-spectral lag relations for seven bursts 
with known redshifts (Norris, Marani & Bonnell 2000) are also plotted for comparison. 
Surprisingly enough, a simple sub-jet model happens to reproduce the observation quite 
well including GRB980425 which has the extremely dim luminosity and large spectral lag. 

The sub-jets model might be compatible with the luminosity-variability relation that 
more complex bursts tend to be brighter (Fenimore & Ramirez-Ruiz 2000; Reichart et 
al. 2000; Schaefer, Deng & Band 2001). According to Plaga (2000), when the power 
density spectrum of the GRB time histories is P oc ~ /^^^^ (Beloborodov, Stern & 
Svensson 1998, 2000; Shaviv & Dar 1995), the variability relates to the pulse width as 
V oc 5'^'^^ oc Wpily^j^. From Figure 2 and §3, the peak flux at a given frequency depends 
on the pulse width as vF^^"-^ oc Wp^jy^f^. Therefore the luminosity- variabihty relation 
is given by ylJp^"-^ oc ~ V^-^iy^'^) for the brightest (dimmest) bursts. In the 

upper panel of Figure 4, we show vLiP^"-^ oc y/(<^+^) by the solid line , using i/LPj^"-^ and 
Wfhwm in Figure 2 and §3 with appropriate off-sets. The observed luminosity-variability 
relations for eight bursts with known redshifts (Fenimore & Ramirez-Ruiz 2000) are also 
plotted for comparison. Surprisingly again, a simple sub-jet model happens to reproduce 
the observation quite well including GRB980425 which has the extremely low variability. 
Although the power density spectrum of the GRBs has been assumed, this agreement 
together with the results of the luminosity-spectral lag relation suggests that a single sub-jet 
model may capture the main features of the GRBs. 

5. DISCUSSIONS 

GRB980425 and its apparent association with SN1998bw first provided a possible 
connection between GRBs and SNe. Since GRB980425 was an unusual GRB, it may be 
a new type of GRB. However, we have shown that the unusually low luminosity, large 
spectral lag and low variabihty of GRB980425 can be interpreted by the off-axis sub-jet 
emission with the viewing angle of 9^ ~ IO/7 ~ 6°. The off-axis sub-jet emission may 
also explain why GRB980425 was not seen in BATSE's highest energy band (> SOOkeV; 
Norris, Marani & Bonnell 2000), since the maximum frequency Pmax becomes lower than 
the observed frequency from Figure 1. Furthermore, the imusually slowly declining X-ray 
afterglow of GRB980425 (oc T^^'^) can also be explained by the off-axis emission from the 
coUimated afterglow with the viewing ang le of ^("^ ~ 30° (Nakamura 1999). Since the axis 
of the afterglow can differ from that of the sub-jet, 9y ~ 6° is compatible with 6*^"^ ~ 30° if 
the sub-jets on the edge of the afterglow were seen. Moreover from the late time light curve 
of SN1998bw, it is suggested that the explosion is jet-like and our hne-of-sight is off-axis 
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from the jet (Maeda et al. 2000). Therefore GRB980425 may be a typical GRB observed 
from the large viewing angle so that our results have strengthened the association of GRBs 
with SNe. 

The bright bursts, such as GRB990123, can be suspected of on-axis emission with no 
or little spectral lag in our model, so that observation with fine time resolution is important 
to measure very small lag. In reality, a burst should have more or less off-axis emission so 
that the lag should exist in the cross-correlation analysis (Norris, Marani & Bonnell 2000; 
Wu & Fenimorc 2000). In this Letter we have studied only the behavior of the single sub-jet 
under the simple treatment of the spectral lag. Nevertheless we have obtained theoretical 
results very similar to the observations, which suggests the viewing angle is one of the most 
important factors to determine the various relations in GRBs. In the next step we should 
consider the multiple sub-jets case using the present results of the single sub-jet.^ 
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Fig. 1. — The total fluence uSi^ in units of 2AorQ7z/Q//5Z}^ is shown as a function of the 
observed frequency by varying the viewing angle < 76*^ < 10, (76*^: integer). We adopt 
^A9 — 1, as — —1, Pb — —2.5 and s = 1. The triangles represent the analytically estimated 
dependence, Umax 6'^ and {uS,^)max 6~^, where 5 — 7(1 — PcosOy) ~ (1 + ^^7^)727 is 
the Doppler factor and the normalization is set by the peak of 7^^ = 10. 
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Fig. 2. — {upper panel): Pulse profiles are shown. The pulse begins at Tgtart = Tq + 
(ro/c/?)(l - /3cos(max[0,^„ - A^])) and ends at T^nd = Tq + (ro/c/?)(l - /?cos(^^ + A^)). 
Here we set Tq = 0. {middle panel): The pulse width WpwHM at half maximum in units of 
r^jcf^'i^ is shown as a function of (1 + 7^6'^). The pulse width Wfwhm can be estimated as 
Wfwhm oc 5 (X {iWOl) for 6^ ~ A^ and Wfwhm oc 5^/^ oc (1+7^^2)1/2 ^ ^ 
which are shown by dashed and dotted lines, respectively, {lower panel): The product of 
the peak flux uPP'^"''' and the pulse width Wfwhm, in units of 2AQrQ'ji'Q/PD'^, is shown as 
a function of (1 +7^^^). This quantity is about the total fiuence uSu in Figure 1. The peak 
flux z/FP""'^ can be estimated from the relation uFp^''''Wfwhm ~ J^S^ oc 
when the maximum frequency Umax is higher (lower) than the observed frequency, which are 
shown by dashed (dotted) lines. In all panels, we adopt 7A^^ = 1, as = —1, Pb = —2.5 and 
s = 1, and the observed frequency is u — 200keV(7i^n/10^keV). 
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Fig. 3. — The isotropic peak luminosity at frequency = 200keV(7Z^Q/10^keV) is shown 
as a function of the spectral lag of the peak time AT by varying the viewing angle 'jOy. 
The corresponding viewing angle j9y is shown in the lower panel. The spectral lag AT is 
defined by the difference of the peak time between frequencies Uj — 200keV(7z/o/10^keV) and 
ux = 20keV(7i/^,/10=^keV). We adopt -fA0 = 1, as = -1, Pb = -2.5, s = 1, ro/cp-f'^ = 10 s 
and an arbitrary vertical unit 2cAo^o7^/-D^. All points with AT < 10~^(ro/c/37^) are plotted 
at AT — 10~^(ro/c/37^). The observed luminosity-spectral lag relations for seven bursts with 
known redshifts are shown from Table 1 in Norris, Marani & Bonnell (2000). Squares and 
triangles arc lags for regions down to 0.5 and 0.1 of the peak luminosity, respectively. The 
dotted lines represent the power law behavior oc AT*^"^"*""^-'/^*"'"^) = AT~^/^ and 

-fOy OC ATV2(^+i) = AT^/i (loka & Nakamura 2001). ^ 
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Fig. 4. — The isotropic peak luminosity at frequency v = 200keV(7Z^Q/10'^keV) is shown as 
a function of the variabihty V by varying the viewing angle 'y6y. The corresponding viewing 
angle jO^ is shown in the lower panel. We adopt jA9 = 1, as = —1, Pb = —2.5, s — 1, 
the power density spectrum of the GRB time histories P oc /'^ ~ f~^^^, and appropriate off- 
sets. The observed luminosity-variability relations for seven bursts with known redshifts 
are shown by filled triangles from Table 1 of Fenimore & Ramirez- Ruiz (2000). The 
analytical estimates of z/L^'^'^*^ oc \/(-2+«B)/(rf+i) _ ^^^peak ^ y(-2+0B)/{d+i) _ 

and z/L^'^"^ oc V^~^^'^'^'^^^^^~^^^ ~ V^^ are shown in each region by dotted, dashed, and long 
dashed lines, respectively. 



